Eigenvalue statistics of reduced density matrix during driving and relaxation 
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We study a subsystem of an isolated one-dimensional correlated metal when it is driven by a 
steady electric field or when it relaxes after driving. We obtain numerically exact reduced density 
matrix p for subsystems which are sufficiently large to give significant eigenvalue statistics and 
spectra of \og{p). We show that both for generic as well as for the integrable model the statistics 
follows the universality of Gaussian unitary and orthogonal ensembles for driven and equilibrium 
systems, respectively. Moreover, the spectra of modestly driven subsystems are well described by 
the Gibbs thermal distribution with the entropy determined by the time-dependent energy only. 
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Introduction. — Spectral universality is one of the key 
features of highly excited complex systems. It has been 
demonstrated and observed in a diverse range of phe- 
nomenologies, ranging from acoustics [H, microwave res- 
onators jif, quantum dots to many-particle systems, 
such as complex nuclei 0] and strongly correlated mod- 
els of condensed matter Q. Universality is quantita- 
tively characterized by the applicability of a parameter 
free random matrix theory (RMT) where the Her- 
mitian operator in question, usually the Hamiltonian, 
is described by an ensemble of Gaussian random matri- 
ces, where the only constraint, whether matrices are real 
symmetric, or complex hermitian, is imposed by the ex- 
istence, or non-existence of a (generalized) time-reversal 
symmetry. Random matrix distribution of energy levels 
is also widely used as a clean indicator of complexity, or 
non-integrability, of a physical model, and is the most 
abstract definition of a quantum chaotic behavior Q ■ 

In this Letter we propose RMT analysis of a completely 
different concept in quantum statistical physics, namely 
of spectra of reduced density matrix (RDM) p of equi- 
librium and non-equilibrium states. We consider RDM 
of strongly correlated quantum systems. In particular 
we study the one-dimensional (ID) model of interact- 
ing spinless fermions (equivalent to a Heisenberg-type 
spin chain) for a variety of simple pure states of the en- 
tire system: the so-called microcanonical (MC) states 
(approximate eigenstates) , the time-evolving states af- 
ter a quench of magnetic flux, or during inductive driv- 
ing with a linearly increasing magnetic flux. We show 
that, quite remarkably, the statistics of eigenvalues of 
RDM of large subsystems is typically described by RMT. 
For equilibrium thermal states, we find agreement with 
Gaussian orthogonal ensemble (GOE), whereas for non- 
equilibrium, driven states, with currents, we find agree- 
ment with Gaussian unitary ensemble (GUE). We note in 
particular, that spectra of RDM of large subsystems typ- 
ically follow RMT even if the entire system is completely 
integrable. 

Furthermore, the RDM can serve as a stringent test 



of thermal properties of nonequilibrium states and their 
thermalization. Our results show that for modestly 
driven systems the entropy density s of the subsystem 
develops in time according to the quasi-equilibrium sce- 
nario, i.e. s depends only on the instantaneous energy 
density e. In agreement with the latter is the observa- 
tion for driven integrable and non-integrable chains that 
the eigenvalue spectra are consistent with the canonical 
Gibbsian form p oc ex-p{—Hcs/T) with well defined ef- 
fective temperature T and HeS being a T— independent 
effective Hamiltonian of the subsystem. 

Model and method. — We study the ID model of in- 
teracting spinless fermions on a chain of even number 
of sites L with periodic boundary conditions. We inves- 
tigate how the system responds to an external electric 
field as introduced in the time-dependent model by the 
varying magnetic flux 4'{t), 

H{t)^-toY.{e^^^'y^ 4ic,+h.c.} 

j 

j 3 

where hj = cjcj, io is the hopping integral, V and W are 
the repulsive potentials between fermions on the nearest- 
neighbor and the next-nearest-neighbor sites. Furtheron 
we use units in which h ^ = to ^ 1. The main 
idea behind introducing W is to break integrability of 
the pure Iq — V model. One expects generic proper- 
ties for the non-integrable case with W ^ 0, whereas 
the integrable system {W = 0) shows anomalous relax- 
ation 8, Jdj and transport characteristics [l2| - [l8| . If not 
stated otherwise, the numerical results for integrable and 
non-integrable cases will refer to V = 1, W = and 
V = 1.4, W ^ 1 systems, respectively, at half-fllling with 
M = L/2 fermions on L = 26 sites. These parameters 
correspond to the metallic regime. 

In the numerical procedure using the microcanonical 
(MC) Lanczos method [l^ we generate initial states 
|5'(0)) for the target energy Eq = (*(0)|ff(0)|*(0)) and 
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the energy uncertainty S'^Eq = (*(0)I[iy(0) -£'o]2|\E'(0)). 
Typically, we consider large Eq (corresponding to high 
T, i.e. /3 = l/T < 0.5) with L = 26 and S ^ 0.01. To 
simulate the MC ensemble, the energy window is small 
on a macroscopic scale (SEp /En <^ 1) but still contains 
a large number of levels |20[. The time evolution of 
^'(0)) — > [9(t)) is calculated by Lanczos propagation 
method [2l| applied to small time intervals {t,t + 5t). 

Since the entire chain is isolated from the surround- 
ings it remains in a pure state In this 
Letter we focus on the reduced dynamics of its subsys- 
tem containing TV subsequent lattice sites. The RDM of 
the subsystem is then p = TrL-Ar|^'(t))(5'(t)| where the 
trace is taken over the remaining L ~ N sites. RDM is 
block-diagonal with respect to number of particles in the 
subsystem n ~ Y^f^^rn ^ Q, N . As the approach is 
numerically accurate the reduced dynamics is exact as 
well. At the same time, the approach allows for sub- 
systems which are sufficiently large to give meaningful 
level statistics as well as spectral and other properties of 
p. With respect to time-dependent response we study 
two kinds of systems: a) driven by a steady electric field 
F = —dt4>{t)/L = const, and b) relaxing but not nec- 
essarily thermalizing after a sudden quench of the flux 
(field pulse): (j){t) = 6'(t)J0. 

Eigenvalues statistics of RDM — We start with the 
presentation of results on the eigenvalue statistics of 
log(p). For the quasi-thermal states p oc exp(— /JiJcff) 
the statistics of eigenvalues of log(p) should be the same 
as the level statistics of the effective Hamiltonian. Since 
we can reach using MC Lanczos method (at half-filling) 
systems with L = 26 the eigenvalue statistics is deter- 
mined for the largest accessible subsystems of = 12 
sites and rt = 6 fermions (containing 924 levels). We 
note that even though the number n of fermions within 
the subsystem is not conserved, RDM is block-diagonal 
with respect to states with fixed n, hence [i?cff, n\ = 0. 
The spectrum {A^} of log(p) is unfolded by a linear inter- 
polation of the integrated density J\f{u>) = J^i ^(w — Ai) in 
intervals containing 40 subsequent eigenvalues. This pro- 
cedure leads to a smoothcned integrated density A/'s,„(a;). 
The unfolded spectrum consists of A; = AfsmiK) where 
we analyze only 2/3 of eigenvalues from middle of this 
spectrum (see [l^l for more details on unfolding). We 
study in the following two standard quantities character- 
izing the level statistics: a) the nearest-level-spacing dis- 
tribution p{x) and the eigenvalue number variance S^((5) 
defined as the variance of the number of unfolded eigen- 
values in the interval of length 5. p{x) and S^(5) char- 
acterize local correlation properties of the spectrum and 
long-range level correlations, respectively. The numerical 
results for log(p) can be compared with the results of the 
RMT for the GOE or GUE ensembles fH,!!!, 




FIG. 1. (Color online) Eigenvalue statistics of log(p) of the 
integrable system: in the equilibrium MC state, in the quasi- 
equilibrium evolution during driving {F = 0.1), and in the 
non-thermal steady state after quenching the flux (Quench). 
Panel (a) shows the integrated spacing distribution I{x) with 
the inset zoom in the low-gap regime, and b) the level number 
variance T,2(5). Lines show results from the RMT. 
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^hoEiS) ^ ^ (^log(2^<5) + 7 + 1 - - ) , 
PGVe{x) ~ exp(-42:^/7r). 



SGUEW^^(log(2^<5)+7+l) 



nx 



exp(-7rx^/4), 



(2) 



where 7 ~ 0.577 is the Euler constant. Note that Hamil- 
tonians of many-body integrable systems have the Pois- 
son distribution with pp{x) = exp(— x) and Sp((5) = 
(5, while generic non-integrable systems with the time- 
reversal symmetry are expected to follow the GOE statis- 
tics. Only cases breaking the time-reversal symmetry 
should result in the GUE statistics. 

We first analyze the spectrum of log(p) for the initial 
MC state in a generic non-integrable system and find 
that both p{x) and (6) accurately reproduce the results 
for GOE (not shown). It is expected since in this case 
without a fiux the time-reversal symmetry is preserved 
and p can be chosen as a real symmetric matrix. In Fig.[T] 
we present numerical data for the integrable system to- 
gether with the prediction of the RMT. The upper panel 
shows integrated spacing distribution I{x) = Ayp{y), 
whereas S^(5) is shown in the lower panel. Surprisingly, 
the eigenvalue statistics of the RDM of a subsystem turns 
out to be independent of the integrability of the total sys- 
tem. 

On the other hand, under a constant (but modest) field 
F > [H Il4l | or after a sudden flux quench |25| we flnd 
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FIG. 2. (Color online) Instantaneous entropy density s{t) vs. 
energy density e{t) for integrable (a,c) and non-integrable 
(b,d) metals. Panels (a,b) show results for systems driven 
by field F = 0.1 and A'^ = 6, 10 subsystems. Dashed curves 
represent high-T analytic results, Eq. (|9]) for L oo. In the 
insets the high-e regime is magnified and 1/L corrections to 
Sao are included and marked as horizontal lines. Panels (c,d) 
show relaxation after flux quenching for A'' = 12 subsystem. 
Here, dashed curves show quasi-equilibrium results for the 
same systems driven with F — 0.1. Arrows mark the direction 
of the processes. 



that the statistics turn into GUE. This is the case for 
the non-intcgrable systems as well as for the integrable 
one, as clearly confirmed in Fig. 1. The GUE statistics 
at F > is consistent with the time-symmetry break- 
ing by a finite current within the subsystem. In the case 
of quenching the decay of the current is not complete, 
at least not within an integrable system where the ab- 
sence of the current relaxation is a hallmark of a finite 
charge stiffness [l2j. We have determined the eigenvalue 
statistics also for the far-from-equilibrium driven states, 
shortly after the electric field has been switched on (not 
shown). The eigenvalues of log(p) are very similar to 
those presented in Fig. [T] for the GUE level case, indicat- 
ing that the RMT statistics of log(p) is very robust. 

Subsystem entropy density- While the von Neumann 
entropy of the total system is at all times zero for the 
pure state \^{t)){^{t)\, the (entanglement) entropy of 
the subsystem S — — Tr7v(plogp) is clearly not. As a 
strong indication of the thermal (or quasi-equilibrium) 
states we can use the relation of the subsystem entropy 
density s = S/N and the energy density e = {H{t))/L 
of the total system (being the same for the subsystem). 
Hence, we present in Fig. [2] the time-evolution of s{t) 
plotted versus £{t) for two systems: driven by a constant 
electric field F = Q.\ (panels a,b) and after a sudden 
flux quench (panels c,d). Results in the right and left 
panels are for the non-integrable and integrable cases, 
respectively. Note first that s{t) only weakly depends on 



N confirming its macroscopic relevance [26|. This is in 
contrast with a specific case of the ground state where we 
have found s(0) oc N^"^ in agreement with the so called 
area laws for the entanglement entropy [23|. 

Since we are studying the high-T regime it is in- 
structive to recall the equilibrium result following from 
straightforward high temperature expansion (HTE) , 

(6) 
(7) 
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where Sco , (^oo refer to T = oo and 1 / L corrections emerge 
due to a restriction of strictly M = L/2 fermions in the 
whole system. Integrating the equilibrium relation /3 = 
ds/de one gets for the equilibrium entropy density 



s(e) = Soo - 



log(2) - 



N - 1 



- 1 



(9) 



(10) 



4L2 6i3 ' 

where again leading 1/L corrections in Soo < log(2) arise 
due to fixing M = L/2. Insets in Figs. [2^, b show that our 
numerical results are very close to this simple estimate 
for Soo- 

The relation (j9l) allows specifying regimes which are 
clearly nonequilibrium or steady but non-thermal. The 
former case occurs, e.g., just after turning on the elec- 
tric field when s = s{£) is convex (see Figs. [2^ and 
[lb) contrary to concave dependence, which according to 
Eq. ([9]) should characterize the quasi-equilibrium evo- 
lution. More interesting is the observation in Fig. [2j: 
that the stationary non-thermal state emerges when in- 
tegrable system relaxes after a sudden quench [1, Q but 
s{t) remains evidently smaller than expected for a ther- 
mal relation, Eq. ([9]). 

Results shown in Fig. [2] may also suggest regimes, 
when the system even out of equilibrium reaches a quasi- 
thermal state. In this case s{e) should become inde- 
pendent of the initial state and close to the prediction 
of HTE. This indeed happens for integrable or non- 
integrable systems driven long enough by a moderate 
steady F (note that F breaks integrability of an inte- 
grable system) or when non-integrable system relaxes 
after the flux quench (see Fig.[2]i). 

Thermal states- More stringent test for a thermal state 
is the requirement that the RDM obeys the canonical dis- 
tribution p oc exp[— /3iJcff (0)] whereby HcS plays the role 
of an effective subsystem Hamiltonian. The effective in- 
verse temperature /3 can be also simply obtained from 
Eq. ([6]). The main open problem concerns the meaning 
of i/oft when the subsystem is strongly coupled to its 
surroundings or is subject to external driving. However, 
we avoid this problem by testing the thermalization hy- 
pothesis without specifying the explicit form of Hcs{(t>). 
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FIG. 3. (Color online) Spectra of log[p(t)] for the integrable 
system driven by field _F = 0.1 determined at time t such 
that e{t) — 0.21 but different initial energy densities eo = 
0.04, 0.1 and different particle subsectors n. Panels (a) and 
(b) shows results for subsystems with A'' = 6 and A'' = 12 
sites, respectively. 
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FIG. 4. (Color online) Spectra of log[p] for the integrable 
system (A'^ = 6, n = 3) driven by = 0.1 determined for 
t corresponding to different /3{t) obtained from HTE. Panels 
(a) and (b) show eigenvalues of log[p] and normalized (and 
shifted) log[p]//3, respectively. 



In order to verify that the quasi-thermal state is de- 
termined solely by the energy density, we have deter- 
mined the reduced dynamics of two identical subsys- 
tems (labelled by subscripts 1 and 2) driven by the same 
field F = 0.1 but starting from different initial energies. 
We compare pi(ti) and ^2(^2) for such times ti and t2 
that both systems have the same instantaneous ener- 
gies (temperatures) E{ti) = E{t2)- Then, one expects 
fog[Pi(^i)] = -PH^s\(t>{h)] + const and log[p2(i2)] = 
— /3i7cff [0(^2)] -I- const. In other words, for the quasi ther- 
mal state operators log[pi(ti)] and log[/92(i2)] should give 
Hamiltonians of the same system but at least subject 
to different fiuxes. Such Hamiltonians may have differ- 
ent eigenfunctions but the energy spectrum should be 
the same (up to eigenvalue fluctuations discussed above). 
Fig. [3] shows results testing the above hypothesis for two 
different subsystems = 6, 12 and two sectors n, respec- 
tively. It is quite evident that the spectra are indepen- 
dent of the initial MC energy density eg, at least for the 
states for which the results on entropy already suggested 
possible thermalization. 

Fig. 13^ finally shows spectra of log[p(t)] for driven in- 
tegrable system (only the largest sector with n = 3 is 
presented). Various curves are obtained for various times 
of driving, when the system has different instantaneous 
energies e[t). Together with Eq. ([6]) e{t) is used to de- 
termine the effective /3(i). For a quasi-thermal evolution 
the spectra of log[p(<)]//?(i) = —i/cff const should 

be the same up to a constant value. As shown in Fig. |4l3, 
even the driven integrable system perfectly fulfills this 
requirement. The thermalization of non-integrable sys- 
tems is commonly expected (apart from a few more spe- 
cific cases [2^, [2^), and our results (not shown) confirm 
that in general. 

Discussion. — Our study shows that the RDM p of 
a subsystem within the isolated system of interacting 
fermions can be a useful tool to investigate proper- 
ties of a nonequilibrium system. Our results on the 



RDM eigenvalue statistics within a ID model of inter- 
acting spinless fermions reveal a universal conclusion, 
that subsystem of an equilibrium MC state obeys the 
GOE eigenvalue statistics, independent of integrability 
or non-integrability of the whole system (note that an 
integrable system as a whole obeys the Poisson statistics 
for the total energy eigenvalues). Moreover, subsystems 
of the driven system and systems quenched with a field 
pulse follow the GUE universality, although the model 
by itself does not break the time-reversal symmetry. 

Further, the spectrum of log(p) contains information 
useful for identifying the quasi-thermal, steady non- 
thermal and the non-equilibrium regimes. For the case 
of quasi-thermal states, which are realized also for finite 
but modest driving, we have demonstrated that the effec- 
tive inverse temperature j3 as the only relevant parameter 
determines the spectra of eigenvalues of \og{p). On one 
hand, this result sets straightforward limits on the relax- 
ation of integrable systems. But possibly more impor- 
tantly, it introduces a nontrivial concept of subsystem's 
effective Hamiltonian TJeff. The physical content and the 
usefulness of the latter, also its relation to the original full 
H has still to be explored. In any case it gives a novel 
approach to investigations of non-equilibrium properties 
of isolated interacting systems, in particular in relation 
to their thermal and transport response. 
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